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ABSTRACT 

The Deep Extragalactic Exploratory Probe (DEEP) is a multi-phase research program dedicated to the study 
of the formation and evolution of galaxies and of large scale structure in the distant Universe. This paper 
describes the hrst five-year phase, denoted DEEPl. A series of ten DEEPl papers will discuss a range of 
scientific topics {e.g., the study of photometric and spectral properties of a general distant galaxy survey, the 
evolution observed in galaxy populations of varied morphologies). The observational basis for these studies is 
the Groth Survey Strip held, a 127 square arcminute region which has been observed with the Hubble Space 
Telescope in both broad I-band and V-band optical biters and with the Low Resolution Imaging Spectrograph 
on the Keck Telescopes. Catalogs of photometric and structural parameters have been constructed for 11,547 
galaxies and stars at magnitudes brighter than 29, and spectroscopy has been conducted for a magnitude-color 
weighted subsample of 818 objects. We evaluate three independent techniques for constructing an imaging 
catalog for the held from the HST data, and discuss the depth and sampling of the resultant catalogs. The 
selection of the spectroscopic subsample is discussed, and we describe the multifaceted approach taken to 
prioritizing objects of interest for a variety of scientihc subprograms. A series of Monte Carlo simulations then 
demonstrates that the spectroscopic subsample can be adequately modeled as a simple function of magnitude 
and color cuts in the imaging catalog. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: distances and redshifts — 
galaxies: evolution — galaxies: structure 


1. INTRODUCTION 

The advent of 8-meter class ground-based optical and in¬ 
frared telescopes, the latest generation of radio and sub- 
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millimeter arrays, and the high spatial resolution of the Hub¬ 
ble Space Telescope (HST) have produced a 100-fold increase 
in observational resources with which to study the evolution 
of galaxies in the local and in the distant Universe. 

The commissioning of the wide-held multihber systems 
used by the 2dF Galaxy Redshift Survey (Colless et al. 2001) 
and the Sloan Digital Sky Survey (Abazajian et al. 2003) 
have hnally enabled the construction of large samples of local 
galaxies, for which a wealth of detailed structural parame¬ 
ters have been determined. Parallel efforts with large optical 
telescopes, leveraged by the 100-fold advantage of the lat¬ 
est generation of multi-object spectrographs, are beginning to 
generate samples of comparable impact in the distant galaxy 
held. These will allow statistically signihcant analyses of the 
evolution of varied populations of galaxies from the present 
to redshifts z ~ 1 (extending over half of the age of the Uni¬ 
verse). 

The Deep Extragalactic Exploratory Probe (DEEP) is a 
multi-phase program focused on the study of the formation 
and evolution of galaxies and of large scale structure across 
this redshift range. The second phase (DEEP2), underway 
at present, uses a sample of ~ 50,000 galaxies with ground- 
based multi-band photometry and spectroscopic redshifts ob¬ 
tained with the recently commissioned Deep-Imaging Mul- 
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tiobject Spectrograph (DEIMOS) on the Keck 2 Telescope. 
This paper describes the first five- year pilot phase (DEEPl). 
The primary observational basis is the Groth Survey Strip 
(GSS) field, a 127 square arcminute region which has been 
observed with HST in both broad I-band and V-band opti¬ 
cal filters and with the Low Resolution Imaging Spectrograph 
(LRIS) on the Keck Telescopes. 

There have been a number of ambitious surveys bridging 
the gap between the local and the distant Universe, and in¬ 
creasing the number counts of faint galaxies over the last ten 
years (for an overview of previous work see Koo & Kron 1992 
and Ellis 1997). Among these works we can cite cf., the ESS 
survey, Arnouts et al. 1997, the LOSS survey, Colless et al. 
1999; the CNOC and CNOC2 surveys, Yee et al. 1996, Yee 
et al. 2000; the Hawaii Deep Eields Survey, Cowie, Songalia, 
Hu, & Cohen 1996. 

Three surveys of comparable size and depth to the DEEPl 
GSS survey are the Canada-Erance Redshift Survey (CERS, 
Lilly et al. 1995), the Caltech Eaint Galaxy Redshift Survey 
(CEGRS, Cohen et al. 2000), and the Very Large Telescope 
Deep Survey with the VIMOS spectrograph (VLTh-VIMOS, 
Le Eevre et al. 2004). Eigure ^ places the DEEPl survey 
in the context similar optical redshift surveys of the distant 
galaxy population, with a direct comparison to the CERS and 
CEGRS data sets. The pioneering CERS is a ^ 600 galaxy 
(250 of which have HST imaging) magnitude-limited survey 
which extends to Iab = 22.5 with a median redshift of z = 0.56 
(Crampton et al. 1995). It is 1.5 magnitudes shallower than 
the nominal limit of DEEP 1, even though the median redshifts 
are similar {i.e., the DEEPl survey focuses on intrinsically 
fainter galaxies), due to differences in the survey sampling al¬ 
gorithms and to common difficulties in obtaining redshifts for 
galaxies beyond z ~ 1.1. The CEGRS survey is focused upon 
the northern Hubble Deep Eield (HDE-N), with a compara¬ 
ble number of redshifts {Zmed = 0.7) for galaxies extending 
to R = 24 in the HDE-N proper and to R = 23 in the Elank- 
ing Eields (Cohen et al. 2000). The recent VLTh-VIMOS sur¬ 
vey centers upon the southern GOODS field (Giavalisco et al. 
2004), containing 784 redshifts for galaxies with HST imag¬ 
ing and an additional 815 in the surrounding areas. The spec¬ 
tral program extends down to I^b = 24, as does the DEEPl 
survey, with a median redshift of 0.73. 

Three features distinguish the DEEP GSS survey in this 
context. Eirst, the field is a continuous band extending over 
40 X 3 arcminutes on the sky (covering a range of 38 co¬ 
moving Mpc by redshift z = 1). Second, the complete HST 
imaging in both V-band and I-band allows for the measure¬ 
ment of a uniform set of rest-frame colors and structural pa¬ 
rameters out to redshifts z ~ 1. Third, the spectral resolution 
of 3-4 A (versus 10 A for the CEGRS and 40 A for the CERS) 
enables both resolution of the [O ll] A3727 doublet feature 
(2.7 A rest-frame split) and that detailed internal kinematic 
measurements be made for individual objects. 

2. OVERVIEW 

We have already published several papers which draw upon 
the DEEPl GSS data set. Initial results of the general red¬ 
shift survey were discussed in Koo et al. (1996), while Vogt 
et al. (1996, 1997) found modest amounts of evolution for 
disk galaxies in the field, Simard et al. (1999) explored the 
effect of surface brightness levels in object detection from the 
HST images, and Im et al. (2001) studied the kinematics of 
massive blue spheroidal galaxies. This paper is the first in 
a formal series dedicated to these DEEP program data. We 



Fig. 1.— The distribution of redshift z versus absolute magnitude, for (a) 
DEEPl (Weiner et al. 2005a), (b) the Caltech Eaint Galaxy Redshift Survey 
(Cohen et al. 2000), and (c) the CERS survey (Lilly et al. 1995; Le Eevre 
et al. 1995; Hammer et al. 1995). The dotted vertical line shows the location 
of M*, and dotted horizontal lines mark z = 1 for each survey. We note that 
the lack of galaxies at redshifts well beyond z = 1 in DEEPl is an artifact 
of limited spectral coverage in the extreme red, rather than being caused by 
a dramatic change in the underlying distribution, while the enhanced counts 
at z = 1 are due to large scale structure across the field {cf., Le Eevre et al. 
1994, Koo et al. 1996). The DEEPl survey goes ~ 1.5 magnitude deeper 
than the CERS, comparable to the depth of Cohen et al., and is distinguished 
by complete HST+WFPC2 coverage. 


outline below the immediate papers published and planned 
within the DEEP Groth Strip Survey Sequence; 

I The Sample (this paper) 

II HST Structural Parameters of Galaxies in the Groth 
Strip (Simard et al. 2002) 

III Redshift Catalog and Properties of Galaxies (Weiner 
et al. 2005a) 

IV Eormation and Evolution of Disk Galaxies from a Sam¬ 
ple of Spatially Extended Velocity Curves (N. P. Vogt 
et al. 2005, in preparation) 

V Evolution of Eield Galaxies in Luminosity and Velocity 
Widths (B. J. Weiner et al. 2005b, in preparation) 

VI Evolution of faint AGN (V. L. Sarajedini et al. 2005, in 
preparation) 

VII The Metallicity of Eield Galaxies at 0.26 < z < 0.82 
and the Evolution of the Luminosity-Metallicity Rela¬ 
tion (Kobulnicky et al. 2003) 

VIII Evolution of Luminous Bulges at High Redshift (Koo 
et al. 2005) 

IX Evolution of the Eundamental Plane of Eield Galaxies 
(Gebhardt et al. 2003) 

X Number Density and Luminosity Eunction of Eield 
E/SO Galaxies at z < 1 (Im et al. 2002) 

Updates to the our mission statement can be found on the 
DEEP website (http://deep.ucolick.orgi, as well 
as a complete data release of the GSS data. 

All Johnson-Cousins magnitudes used in this paper are on 
the Vega system (with the exception of those shown in Eig- 
ure|3 which uses the AB system); transformations between 
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Fig. 2.— The distribution of objects along the GSS. The top two panels show the extent of the strip on the sky, running from the first chevron (number 4, at 
the northeastern extreme) to the 28'* chevron (number 31) and covering an area roughly forty by three arcminutes. The upper panel shows the distribution of all 
objects detected in either or both and V '606 passbands in the FOCAS catalog, and the lower panel the subset for which we have obtained Keck+LRIS spectra. 
The lower panels show the fraction of targets within the FOCAS catalog down to the magnitude limit of our general redshift survey for which we have spectra, 
as a function of position along (left panel) and across (right panel) the GSS. 


the Vega and AB systems are listed in Appendix B. We have 
chosen to use a ACDM cosmogony with h = 0.1 = 0.3, 

and = 0.7 for this entire series of papers, unless otherwise 
noted. 

3. FIELD HISTORY 

The GSS held is a high latitude held in the HST continuous 
viewing zone with no bright stars and with marginal levels of 
galactic extinction. The primary science drivers for the initial 
HST+WFPC2 Cycle 4 observations were to study cosmology 
and galaxy evolution by evaluating the evolution of number 
counts, the luminosity function, colors and morphology, and 
the distribution of galaxies and of clusters in redshift. 

The held is dehned as the spatial position of a long 
band of imaging data from two Hubble Space Telescope 
and Wide Field Planetary Camera (HST+WFPC2) programs 
(GTO 5090 with PI E. Groth; GTO 5109 with PI J. West- 
phal). They consist of 28 overlapping WFPC2 chevrons, num¬ 
bered 4 through 31 because three chevrons originally placed 
at the northeast end of the strip were discovered to lie too far 
from appropriate guide stars and were replaced with three new 
pointings at the south-west end. The held is oriented at 45° 
eastwards of north, centered at (14*17"',-1-52°), as shown in 
Figure|2] has a total area of 127 square arcminutes, and lies at 
a Galactic latitude of h ^ 60°. 

Though the strip of WFPC2 chevrons runs along an angle 
of 45° eastwards of north, the orientation of the roll angle, 
equivalent to the axis of symmetry of each chevron along the 
line bisecting each chip 3 from pixel (0,0) to (800,800), lies at 
40.5° eastwards of north. Because of the 4.5° offset, chips 2 
of each chevron overlap slightly within a region of size 19" by 
21" with the chips 3 of the next chevron along the strip. These 


overlap regions were used to align the chevrons together in to 
a single supermosaic for the data taken within each passband, 
as were diffraction spikes that crossed chip to chip bound¬ 
aries. 

Each chevron has been imaged for 2,800 seconds in a broad 
V hlter (E606W, Veoe) and for 4,400 seconds in a broad / 
hlter (E814W, /gu), with the exception of the single deep 
chevron (number 7, fourth from the northeast edge, centered 
at 14*17'"30‘', -f 52°30"') which was observed for 24,400 sec¬ 
onds in Veoe and for 25,200 seconds in Igu- Images for the 
primary program were constructed by rejecting cosmic rays 
from four undithered exposures in each filter, each lasting 700 
seconds in V or 1100 seconds in /. (A series of four orbits 
were spent observing each chevron, with each orbit consist¬ 
ing of a single exposure in V and a single exposure in I.) Eor 
the deeply imaged chevron, eight 2100 second and four 1900 
second exposures, each taking a full orbit, were combined in 
y, and twelve 2100 second exposures were combined in I. 

Additional WFPC2 optical imaging data have been col¬ 
lected along the GSS by a variety of authors and teams, in¬ 
cluding 3 X 2,800 seconds in a broad B hlter (E450W, B 450 ) 
and 3 x 2,800 seconds in /gi 4 (GO 5449 with PI S. Lilly) for 
a chevron which overlaps by 70% with the third GSS chevron 
(number 6 ). In Cycle 9, all GSS chevrons but the deeply im¬ 
aged number 7 were re-imaged for 700 seconds with Veoe (GO 
8698 with PI J. Mould). A small number of objects within 
the Groth Strip have been selected for additional observations 
with other HST instruments (NICMOS GO 7871 with PI A. 
Connolly; NICMOS GO 7883, STIS GO 10249 with PI N. 
Vogt). 

Keckn-LRIS ground-based BRI imaging of the northeast 
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one-third portion of the Groth Strip has also been conducted 
by the DEEP team. Shallower UBRl imaging (down to U = 
25, B = 25.5, R = 23.4,1 = 23.8 for S/N = 10) of the same re¬ 
gion was done at KPNO, and may be extended as a part of the 
DEEP2 survey. Additional data have been obtained at longer 
wavelengths for parts of the field using SCUBA and the VLA 
(Eomalont et al. 1991), as well fairly extensive data in K-band 
(c/, Cristobal-Hornillos et al. 2003; Conselice et al. 2005). 
The field also has been or will be observed as a part of numer¬ 
ous large survey programs (e.g., the CEHT Legacy Survey; 
also with XMM, Chandra, Galex, JWST, and ISO). This list 
is not a complete census; we simply wish to illustrate that a 
significant amount of observational resources, spanning mul¬ 
tiple wavelength regimes, have been devoted to this portion of 
the sky. 

4. ASTROMETRIC SOLUTION 

An initial astrometric and photometric analysis of the GSS 
WEPC2 data was completed with the Eaint Object Classifica¬ 
tion and Analysis System (EOCAS; Jarvis & Tyson, 1981) for 
11,547 objects down to 33''“^ magnitude in (Veoe + Jsm)/2. 

An independent astrometric solution was created for a su¬ 
permosaic in each passband composed of all 28 chevrons, by 
first transforming the individual chip coordinates to a com¬ 
mon coordinate system based on the polynomial coordinates 
derived by Holtzman et al. 1994. The rms residual errors in 
X and in y on each chip are equal or less than 0.02" along 
the entire strip, as determined by evaluation of the objects ob¬ 
served on two chevrons in the overlap regions. The values are 
fairly uniform, though we note that the residuals on the chip 
3 regions are slightly higher on average (by 25%) than those 
for chips 2 and 4. 

There are six HST guide star catalog objects within the 
strip, which were used for astrometric zero pointing. The er¬ 
ror in guide star coordinates is a sizable fraction of an arcsec- 
ond (we estimate 0.5" in each coordinate), which will affect 
the absolute calibration. Star positions were measured on the 
WEPC2 chevrons by measuring the intersections of diffrac¬ 
tion spikes; in one case the stellar flux peak could also be 
measured, and the two techniques agreed to within a pixel. 
An rms residual of 0.33" was obtained by using all six stars, 
and the best, final fit was obtained from four stars (rms value 
of 0.11"), by discarding two stars which fell very close to the 
edge of the WEPC2 chips. 

Because of their late attachment to the Cycle 4 observa¬ 
tional program, we note that chevrons 28-31 were taken with 
an offset in roll angle of 0.6004° relative to the remainder of 
the strip (the roll angle was specified to only 0.01° in the orig¬ 
inal observing proposal). In contrast, the average error in the 
roll angle was measured to be 0.2° across the entire Cycle 4 
WEPC2 data set. These factors were taken into account in the 
astrometric solutions. A comparison of the final solutions for 
the Veoe and/gn data yielded errors of less than 0.01" between 
the two passbands; this close agreement allowed us to overlay 
the two mosaics on top of each other to produce a valid color 
map of the entire strip. 

Summarizing the primary sources of error in the astrometric 
model, small scale (i.e., chevron scale) errors are controlled 
by the chip distortion model (0.02" rms), chevron offset de¬ 
termination (0.025" rms), and roll errors (0.025" rms in the 
chip corners furthest from the chevron center), leading to a 
total rms error of 0.035". This translates into rms errors of 
order 0.1" in right ascension and declination, and a rotational 
uncertainty of 0.017°, along the entire strip. When evaluating 


individual objects, one must of course also take into account 
the errors in object centroiding (discussed below). 

5. PHOTOMETRIC CATALOGS 

The photometric catalog was created using EOCAS within 
the Image Reduction and Analysis (IRAE) environment to 
determine valid detections, classified as 10,955 galaxies, 436 
stars, and 156 objects of uncertain morphology (e.g., fuzzy, 
diffuse) where classification was done by rescaling (broaden¬ 
ing and narrowing) the point spread function (PSE) to create 
a set of templates then used to best match each object’s flux 
distribution. The objects identified as stars span the full range 
of magnitude of the sample, and represent a conservative esti¬ 
mate of the stellar contribution. The spectral subsample con¬ 
tains 37 objects classified as stars; divided between 28 objects 
at z = 0 and nine for which a redshift could not be determined 
(due to e.g., faint flux, lack of distinguishing features). An 
additional 128 spectra were obtained for stars which had been 
misidentified as galaxies by EOCAS. In summary, EOCAS 
successfully identified from 55 to 73% of all stars brighter 
than 19'* magnitude, dropping to the range of 10 to 22% for 
all magnitudes, within the spectroscopic subsample. There 
was in general great difficulty in distinguishing between in¬ 
termediate redshift (z ~ 0.7) compact galaxies and stars, e.g., 
halo M dwarfs, solely from photometric data. Note that the 
EOCAS object classification codes were not used during the 
process of selecting spectroscopic subsamples. 

EOCAS operates by detecting objects of a set minimum size 
that register above a set threshold level above a sky back¬ 
ground level. Background levels were determined for each 
chip (2 through 4, note that we did not use the planetary cam¬ 
era chip data, so as to maintain a uniform resolution and sen¬ 
sitivity throughout) and on each chevron, by examining blank 
rectangles of sky. The detection threshold was set to 3a above 
sky background for the bulk of the GSS, and to 4.5cr for the 
deeply imaged field (chevron 7). When making this evalu¬ 
ation, the minimum object size was set to two pixels (0.2") 
across the entire strip. 

Three objects were manually added to the catalog of EO¬ 
CAS detections, each being a target that fell well within the 
extended flux of a brighter galaxy or stellar diffraction spike 
and thus was not recognized as a separate object, but for 
which we wished to obtain a redshift nonetheless. 

After each object detection, EOCAS computes an individ¬ 
ual sky background surrounding each object. Each object was 
then “grown” by adding concentric rings of pixels until its to¬ 
tal area had doubled. The grown object was then evaluated 
to determine its magnitude and centroid. EOCAS attempts to 
split blended objects automatically, but the algorithm used for 
this purpose was not found to be optimized for faint object 
detections on these images (i.e., down to (V 606 + 78i4)/2 = 25 
completeness limit, as discussed below). It was thus neces¬ 
sary to examine each object by eye, to verify and on occasion 
to correct the division of parent objects in order to produce 
the proper splitting into de-blended objects. The low detec¬ 
tion threshold and small object size also introduced spurious 
“noise” objects, which were filtered out by matching Veoe and 
/814 frame detections. 

One potential drawback to this technique is that the size of 
the detection area for a given object can vary between Veoe and 

** IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in Astron¬ 
omy (AURA) under cooperative agreement with the National Science Foun¬ 
dation. 
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TABLE 1 

Catalog of Aperture Photometry 


Identification Codes 


GSS 


B" 


CFRS" 


Coordinates (J2000) 
R.A. Dec. 


V. 


Magnitudes 

Vc-Ic, Vf If 


V, 


Selection 
Ig Criteria 


072_4040 G82a B08839 
073_0542 G1307 B09464 
073_1810 G1318 B09321 
073 2356 G1355 B09831 


14:17:37.836 +52:28:58.45 22.15 20.92 1.27 20.64 19.63 19.58 20.61 disk 

14:17:46.405 +52:28:39.65 24.70 23.63 1.08 24.59 23.55 22.06 23.43 zsurvey 

14:17:42.635 +52:28:45.32 23.48 22.01 1.43 21.96 20.54 21.30 22.85 morph 

14:17:47.401 +52:29:00.50 23.51 22.89 0.64 23.42 22.82 22.55 23.09 serendip 


073_3539 G1367 B09751 C14.1043 14:17:45.385 +52:29:08.27 22.04 20.35 1.77 21.19 19.65 19.60 21.12 phz 


"See Rhodes, Refregier, & Groth (2000); '’see Brunner, Connolly, & Szalay (1999); "see Lilly et al. (1995); ‘'HST F606W 
andF814W magnitudes, Veoe and/gw within l."5 diameter, 1406-/814 within l" diameter; Veoe and/gu measured with 
FOCAS; V4o6 and /gi4 measured with GIM2D. 


/814 frames, as the data for the two bandpasses have different 
levels of sensitivity. In addition, an actual variation in object 
size can be a function of Veae-hu color, one of the key cat¬ 
alog parameters which was used in selecting a subsample for 
follow-up spectroscopic observations. Object centroiding was 
also done independently in each bandpass, leading to shifts in 
object position between the two bandpasses. Taken together, 
these factors might strongly affect the applicability of the “to¬ 
tal” magnitudes derived from FOCAS for our fainter candi¬ 
dates. 

We note that the Veoe catalog data appear to extend more 
deeply than those of the /814 observations {i.e., more objects 
are detected in Vgoe but not in /g^ than the reverse). This effect 
becomes significant for objects fainter than 25 in Vgoe, beyond 
the limits of the spectral sample. Down to 25 in Veoe and to 24 
in /814 (i.e., extrapolating a median Veoe-hu color of one for 
non-detections), there are 171 objects detected only in Veoe 
and 157 more detected only in /g^ (and within the spectral 
sample only 16 targets were undetected in /gi 4 and 23 were 
undetected in Veoe)- In contrast, when we extend these limits 
by one magnitude we find 1681 objects detected only in Veoe 
versus 335 detected only in /g^, and the numbers increase to 
5198 versus 760 when the entire optical catalog is evaluated. 

This increase in faint Veoe detections is caused by several 
factors. First, object detection was set to penetrate well into 
the noise, and the faint detection threshold was set indepen¬ 
dently in each bandpass and extends to fainter levels in Veoe- 
Second, the narrower Veoe point spread function (cf., Caser- 
tano et al. 2000) allowed “blobby” objects to be split into 
multiple components, while the same group of objects (or, 
alternatively, the same galaxy plus H II regions) was treated 
as a single object with the summed, brighter magnitude in the 
/gi 4 data. 

We defined a set of object names according to the following 
convention. Every object within the imaging catalog received 
a seven-digit name, optionally ended with a single letter. The 
first three digits are the WFPC2 chevron number along the 
strip (two digits between 4 and 31, with a leading zero if nec¬ 
essary) and chip (single digit, one of chips 2 through 4); they 
are followed by an underscore symbol. The final four digits 
are the x and y position of the object on its chip as detected 
in /gi 4 , if so detected (else in Veoe), rounded off to the nearest 
ten pixels. For example, an object found in field 7 (the deeply 
imaged field) on chip 4 at coordinates (235, 516) would thus 
be named 074_2452. 

One can conceivably (though rarely) have multiple objects 


with the same name (when objects are found within the same 
1" square on a chip), so we also track the FOCAS-generated 
identification number (running from 1 to 11,543) to avoid any 
confusion. When we do include two such objects in the spec¬ 
tral sample, the second to be observed has an "a" attached to 
the end of its seven-digit name. 

We elected to conduct a second photometric analysis of the 
data based on aperture magnitudes, measured with the IRAF 
package APPHOT. Object centroids as determined from the 
FOCAS /814 data were used to define the location of all tar¬ 
gets within both bandpasses, supplemented by the Veoe object 
centroids for those undetected in /g^. Aperture magnitudes 
were determined for both 1.0" and 1.5" diameter regions for 
each object detected by FOCAS. 

We analyzed 421 objects brighter than 25'* magnitude 
within the overlap regions of successive chips 2 and 3 twice 
(i.e., once per separately observed chip, or chevron) to deter¬ 
mine rms error rates. Object centroids had been fit to the chip 
3 data with the FOCAS catalog; these were used to determine 
object positions on the overlapping chip 2 regions. The offset 
between successive chevrons varied by ±2" across the en¬ 
tire GSS, and was determined for each chevron by examining 
the most compact objects in the overlap region. No rotation 
correction was applied, as the linear offsets were deemed suf¬ 
ficient. The individual derived chip 2 object positions were 
then allowed to shift to match the local maximum, resulting 
in sub-pixel corrections. Background levels were estimated 
on a chip-by-chip basis, and then aperture photometry was 
performed using set 1.0" and 1.5" diameter apertures for both 
sets of images. 

The differences in flux levels were found to be fairly con¬ 
stant, at 0.025 ± 0.030 magnitudes throughout the /g^ data 
set. Comparison values for the 1.0" and 1.5" apertures dif¬ 
fered by 0.003 magnitudes, and there was no significant trend 
with magnitude for objects brighter than 25'* magnitude. This 
is likely a combination of two competing effects. The brighter 
objects exhibit high flux levels but extend spatially beyond the 
aperture limits, and so small variations in centroiding between 
chips 2 and 3 can cause variations in estimating the aperture 
flux. The fainter objects are contained almost entirely within 
the apertures, and thus centroiding issues become less rele¬ 
vant, but the overall signal to noise levels decrease. The Veoe 
data showed slightly larger differences between chip 2 and 
chip 3 measurements, at 0.035 ± 0.035 magnitudes, due in 
part to slightly elevated background levels. 

For spectroscopic target selection, we elected to use the Veoe 
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and /814 magnitudes as measured in the 1.5" diameter aper¬ 
tures, combined with the 1.0" diameter aperture Veoe- /si 4 
colors (“Core” colors) to define object magnitudes and col¬ 
ors. This improved the accuracy of color measurements (i.e., 
a relative flux measurement), particularly for faint detections, 
and permitted us to extend color and two-bandpass measure¬ 
ments to faint objects that had been detected by FOCAS in 
only one of the two HST bandpasses. We label this catalog 
of magnitudes and colors the “Core” catalog, as the object 
fluxes are sampled only in the core regions {i.e., aperture size 
is constant, and does not vary with object angular size). 

Table ^ shows a representative set of entries to the ob¬ 
ject catalog, available in its entirety on the DEEP web site 
at http: / /deep. ucolick . org Each object is identi¬ 
fied with its DEEPl seven to eight element name, followed 
by the EOCAS catalog identification and the names (as avail¬ 
able) in the ground-based UBRI (Brunner, Connolly, & Sza- 
lay 1999) and in the CERS survey object catalogs. Coordi¬ 
nates are listed (J2000), followed by aperture Vgoe and /g^ 
aperture magnitudes and Veoe” ^814 colors. The final column 
indicates the selection criteria by which the object was placed 
on a Keck-tLRIS mask, if this occurred (see discussion of key¬ 
words below). 

A third photometric analysis was conducted with the galaxy 
modeling package GIM2D after the majority of the spectral 
program had been completed, to derive a detailed, extended 
set of structural and morphological parameters for all objects. 
Discussion of this “GIM2D” catalog, and its optimization, is 
the primary topic of Paper II (Simard et al. 2002). These 
data, including magnitudes as well as structural parameters, 
are used throughout most of the DEEPl science papers. 

Pigures|3and|4]demonstrate the relation between photomet¬ 
ric measurements made through these three analyses. We em¬ 
phasize that, though the photometric analyses are independent 
for each catalog, the object centroids were all drawn from the 
EOCAS catalog and were not relocated in subsequent analy¬ 
ses {i.e., we have preserved the link between an given object 
name and its position of the WEPC2 chevrons, within all three 
catalogs). The “Core” magnitudes are seen to systematically 
underestimate the fluxes of the brightest objects, as expected 
given that their angular sizes extend well beyond the 1.5" di¬ 
ameter aperture. Note that the effect is most significant for 
objects brighter than 23'’'^ magnitude, which were all placed 
in a single magnitude bin during our spectral target selection 
process, and does not vary significantly with V-I color. When 
comparing the two independent estimates for total magnitude, 
we find that the dominant source of variation is the scatter in 
the measurements at faint (below 23'^'^) magnitudes. 

There are a few {N ~ 20) objects with large offsets from 
catalog to catalog; these extreme differences, on the level of 
two to five magnitudes, are primarily caused, in order of fre¬ 
quency and importance, by (a) difficulties in measuring fluxes 
for objects on the extreme edges of chips (centered within 
the outer 7", where the background measurement has been 
compromised by the sharp decline in pixel responsivity and 
differences in estimating background levels assume a criti¬ 
cal role), (b) overlapping objects which are not deconvolved 
identically by each set of algorithms, and (c) stellar diffrac¬ 
tion spikes. We examined all such outliers, and found only 
one which could not be explained by one of these three fac¬ 
tors: 063_2764, an extremely diffuse, extended low surface 
brightness object. 

The FOCAS and Core catalogs show the fewest such off¬ 
sets, as shown in the first three panels of Figure 0 The 



Fig. 3.— The dispersion of magnitudes throughout the three imaging anal¬ 
yses (FOCAS total, 1.5" aperture Core, and GIM2D total) used for the GSS 
sample. The first three panels show the relative distribution of magnitudes be¬ 
tween catalogs; the dashed line follows the relation y = x and the solid line is 
an unweighted fit to objects between 19 and 25 magnitudes (boxed) differing 
by less than 2 magnitudes between catalogs. The final panel plots the dif¬ 
ference between FOCAS and Core magnitudes versus the difference between 
GIM2D and Core magnitudes, with a fit to the same points as in the other 
panels. We find that the Core magnitudes systematically underestimate lumi¬ 
nosity for the brightest objects, because their large angular diameters extend 
beyond the 1.5" Core apertures, while the difference between FOCAS and 
GIM2D total magnitudes becomes more significant for objects fainter than 
23'^^ magnitude. The final panel demonstrates that the offset for the Core 
aperture measurements is more significant than the scatter between the two 
estimates based on total magnitudes, as the data tend to be forcibly extended 
along the y = x line by the Core offsets (with the brightest objects extending 
this trend preferentially to the lower left hand quadrant) rather than scattering 
uniformly about the origin. 


GIM2D catalog is more sensitive to the first, and most com¬ 
mon, cause of this type of error because its algorithms uti¬ 
lize a larger region around each object when determining the 
background flux level. It is thus likely that the GIM2D magni¬ 
tudes are the least accurate measure of flux for an object when 
a substantial (greater than two magnitudes) difference exists 
between catalog entries. 

When examining relative Veoe and /giq magnitudes {i.e., col¬ 
ors), we found that the independent centroiding between pass- 
bands in the FOCAS measurements produced a significant 
scatter in color overall, and an increase in colors at the ex¬ 
tremes of the color distribution (particularly on the blue end). 
We describe two representative cases for clarity. In the first 
case FOCAS identified two objects separated by 0.9" on the 
sky in/ 814 , one of magnitude 21.8 and the other 25.3. In Vgoe, a 
single object with magnitude 22.6 was placed within the same 
region. The single bright Veoe detection was then matched to 
the fainter /giq object, while the brighter /gu object was as¬ 
sumed to have been undetected in Veoe- The resultant mis¬ 
match produced an object with an artificially low /gi 4 flux, 
resulting in an extremely blue color index. In the second case 
FOCAS split a “blobby” object into three separate objects on 
the Veoe image, while melding the total flux into a single ob¬ 
ject on the /814 frame. The match of the total object in /gi 4 
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Fig. 4.— The distribution of colors throughout the three imaging analy¬ 
ses (FOCAS total, GIM2D total, and 1.0^^ aperture Core) used for the GSS 
sample. The first three panels show the relative distribution of V 606 “ 4i4 
color between catalogs; the dashed line traces the relation y = x and the solid 
line is an unweighted fit to objects between 19 and 25 magnitudes differing 
by less than 2 magnitudes between catalogs (no significant change appears 
in the fits if we relax these magnitude limits). The final panel plots the dif¬ 
ference between GIM2D and FOCAS colors versus the difference between 
Core and FOCAS colors, with a fit to the same points as in the other panels. 
The FOCAS-derived colors contain some substantial outliers, particularly in 
a population of objects with extremely blue measurements of color which are 
not reproduced in the other two catalogs. The final panel demonstrates that 
this offset is more significant than the difference between the GIM2D and 
Core estimates of colors (which are offset by less than half a magnitude), 
as the data tend to be forcibly extended along the y = x line by the FOCAS 
offsets rather than scattering uniformly about the origin. 


versus a single component in Veoe produced an extremely red 
color index. These outlying cases can be characterized by to¬ 
tal object areas which differ greatly (by factors of 10 to 90) 
between the two frames (versus a factor of ^ three across the 
rest of the data set). These objects are mainly found in the lu¬ 
minosity range between 23'’'^ and 25'^' magnitudes. Extremely 
blue and extremely red objects were rated to be of consider¬ 
able importance within our spectroscopic program, thus it was 
a priority to determine accurate colors for them. 

The use of a single set of object positions for both the Veoe 
and the /gu images in the Core and GIM2D catalogs pro¬ 
vides a more robust estimation of the relative flux (though 
the background levels were still set independently within each 
passband). We note that the 1.0" aperture Core colors are 
0.03 ± 0.15 magnitudes redder than those determined for the 
entire galaxy with the GIM2D package, where the positive 
offset could be due to a slight nuclear concentration of redder 
stellar populations in many objects. 

One of the science projects driving the spectral program 
was the modeling of photometric redshift estimates based on 
multi-color broad-band photometry (Brunner 1997; Brunner, 
Connolly, & Szalay 1999). We were actively searching for 
a small population of faint galaxy candidates with extremely 
blue colors, to use to improve the models, as well as trying 
to study the evolution of extremely blue and red objects. The 


color information used for sample selection was thus deemed 
to be quite important. For this reason, we elected to use aper¬ 
ture colors, and magnitudes, when creating the spectral sam- 





Bulge Fraction Bulge Fraction 


Fig. 5.— The distribution of number counts and bulge fractions as a func¬ 
tion of hi^AB magnitude. Dotted lines on the top panels show the three 
bins used in the spectroscopic target selection, with a fourth bin for fainter 
objects extending down to = 26. The top left hand panel shows the 

distribution of number counts for the 9212 objects in the GSS imaging cata¬ 
log and the 522 within the HDF-N field imaging catalog (Marleau & Simai'd 
1998). The HDF-N catalog extends considerably deeper at faint magnitudes 
{huAB 25) than the CjSS catalog, due to the extiemely deep, repeated 
exposures of the field. However, there is quite good agreement in the dis¬ 
tribution within the limits of our spectral sample (Veoe+ f 8 i 4)/2 ^ 25 (note 
that (V6oe+f8i4)/2 in the Vega system is roughly equivalent to /gu in the AB 
system). The top right hand panel compares the number counts for the GSS 
spectral sample with the HDF-N field catalog. There is reasonable agreement 
in the distribution within the limits of the brightest bin of our spectral sample 
where (V 606 + 4i4)/2 < 23. The bottom panels compare the bulge fraction 
distribution within the HDF-N imaging catalog and the GSS spectral sample, 
finding an excess of no bulge component objects in the HDF-N sample which 
is concentrated at faint magnitudes and does not hold for objects brighter than 
huAB = 2^. 


Figure 13 evaluates the depth of the HSTh-WFPC 2 imag¬ 
ing, and the consequent limits to the photometric catalogs, 
by comparing two analyses of /g^ data conducted with the 
GIM2D package. We compare the distribution of number 
counts (Marleau & Simard 1998) within the drizzled and 
deeply imaged HDF-N field (Williams et al.l996) with those 
within the GSS imaging catalog, scaling the absolute number 
of counts (522 versus 9212, for objects brighter than /gi 4 = 26) 
by the relative surface areas (4.5 versus 127 square arcmin- 
utes), and find good agreement down to /814 = 25, with major 
differences becoming apparent beyond 25 magnitudes. Note 
that Veoe data were not used in the GIM2D HDF-N analysis 
and are thus not available for comparison, and also that the 
HDF-N I-band magnitudes were cataloged under the AB sys¬ 
tem. A parallel comparison with the 785 objects in the GSS 
spectral sample with GIM2D parameters yields more scatter, 
but the agreement is still reasonable for objects brighter than 
(V606 + f814)/2 < 23. 

The distribution of bulge fractions within the HDF-N sam¬ 
ple contains an increase in objects best fit with no bulge com- 

















































8 


ponent, relative to both the GSS imaging and the GSS spec¬ 
tral samples. This effect is strongest for faint objects, and is 
not significant below 23 magnitudes. This is both because 
the difference in bulge fraction distributions is not present in 
the brighter objects, and because the number of objects has 
dropped to 72 for HDF-N sample, decreasing the statistical 
significance of the comparison process. Given the catalog 
limits, we speculate that these objects may be low surface 
brightness late-type spiral galaxies, falling below the surface 
brightness sensitivity of the GSS samples. 

In summary, we estimate the 5cr detection limit of the imag¬ 
ing catalogs to lie at (Veoe+78i4)/2 = 25 (though note that 
surface brightness completion becomes an issue for signifi¬ 
cantly brighter objects, at roughly 23 magnitudes). Note that 
errors on individual magnitudes are 0.02 for observed quan¬ 
tities, and 0.02 - 0.03 for derived rest-frame magnitudes (see 
Simard et al. 2002 for complete details of such). 

6 . SAMPLE SELECTION 

A series of spectroscopic samples was then selected, using 
the Core (1.5" diameter aperture) magnitudes and the Core 
(1.0" diameter aperture) colors to define the distribution of 
fluxes. Targets were selected to populate three magnitude 
bins, and divided into five color bins as well, for a general 
redshift survey. Table |2] shows the distribution of the com¬ 
pleteness ratio (fraction of objects from the imaging catalog 
within each bin targeted for spectroscopy) for the total spec¬ 
tral survey program, which includes a set of objects chosen 
for additional properties (described in detail below) as well as 
the general redshift survey. 

The primary magnitude limit for the general survey was 
(V 606 + 7814)/2 = 24, well above the detection limit of 25 for 
photometric completeness (though see Simard et al. 1999 for 
a discussion of the effect of surface brightness biases). This 
limit was chosen to lie well above the completeness level of 
the optical catalogs, as established by comparison with the 
data for the deeply imaged HDF-N field (see Figure|3i, and to 
obtain sufficient spectral flux for redshift determination of ob¬ 
jects exhibiting a range of morphologies (e.g., both absorption 
and emission dominated spectra). 

There are 83 objects with magnitudes between 24 and 25 for 
which spectra were obtained, for the following reasons. The 
pool of viable targets was deliberately widened to 25 magni¬ 
tudes when sampling chevron 7, to extend the spectral survey 
to fainter magnitudes in the region where our photometry was 
significantly deeper than the norm. This resulted in observa¬ 
tions of 28 objects which lay below the primary magnitude 
limit. Six additional faint objects were included as they were 
high priority targets for special subsamples (R#, hi-z, or phz, 
defined in Table |3- Twelve more “filler” objects were in¬ 
cluded due to a lack of brighter previously undetected spec¬ 
tral candidates within a given “slit position" along the GSS, 
serving as candidates of last resort for the north-eastern end 
of the GSS where our spectral observations were most con¬ 
centrated. Finally, 37 objects with magnitudes between 24 
and 25 were obtained serendipitously, falling within slits de¬ 
signed for brighter targets. There are also six fainter serendip¬ 
itous objects within the spectral sample, lying between 25'* 
and 25.5'* magnitude. 

Note that Table |2] contains data on 788 objects from the 
spectroscopic sample, omitting 30 objects for which Core 
magnitudes could not be measured due to object close prox¬ 
imity to a WFPC2 chip edge or interfering stellar diffraction 
spikes. In addition, five stars selected for use in the mask 


alignment process are contained within this sample that were 
not listed in the spectroscopic sample described in Weiner 
et al. 2005a. 

Figure|6lshows the distribution of the imaging and the spec¬ 
tral samples in color and in magnitude. We prioritized ob¬ 
jects with blue colors (Vgoe-fsM < 0.5) or with red colors 
(V606“7814 > 1 -75) in the selection process; they are thus over¬ 
represented in the spectral sample at a higher fraction (150%) 
than objects with intermediate colors in the range above 24'* 
magnitude. Because the five color bins were chosen to be 
symmetric about the total color distribution, however, this did 
not significantly affect the relationship of average color as a 
function of magnitude in the spectral subsample. In both sam¬ 
ples, average Vgoe - 7814 colors are relatively constant for mag¬ 
nitudes brighter than 23 and then begin to shift blueward, at a 
rate of 0.2 magnitudes in color per unit of total magnitude. 




Fig. 6 .— The distribution of Core color (1.0" diameter aperture) versus 
Core magnitude (1.5" diameter aperture) for the magnitudes used in spectral 
target selection, for (a) the complete imaging catalog and (b) the spectral 
subsample selected to observe with Keck+LRIS. Colors have been averaged 
in half-magnitude running bins, and the average values are connected with 
a line against the field of objects to show the mean trend. The dotted lines 
divide the space into fifteen blocks brighter than (Im + Vf/ad)/'! = 25; each 
bin was tracked separately through the spectral target selection process. 


7. SPECIAL SUBSAMPLES 

One could say that the spectroscopic program was designed 
in an umbrella factory, in the sense that a host of scientific sub¬ 
programs motivated the inclusion of various subsamples of 
objects, chosen with additional criteria beyond the described 
optical magnitudes and colors. In addition to the general mag¬ 
nitude limited sample defined above, several types of objects 
were prioritized in small numbers for selection from the pho¬ 
tometric catalog. Table 0 gives a brief description of these 
additional selection criteria and their fractional representation 
within the complete spectral sample. We emphasize that al¬ 
most all of these objects were also eligible for selection as a 
part of the general redshift survey (the exception being six 
serendipitously observed extremely faint objects which fall 
beyond the magnitude limit). 
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TABLE 2 

Observed Sample Completeness 




Color 

m < 23 

m 

23 < m < 24 

= {V+I)/2 

24 < m < 25“ 

25 < OT < 25.5* 



V-I < 

0.25 (bluest) 

1/2 

0/2 

1/99 

1/194 

0.25 

< 

V-I < 

0.50 (blue) 

6/16 

29/90 

10/454 

imi 

0.50 

< 

V-I < 

1.75 

333/941 

229/1354 

60/2613 

3/1372 

1.75 

< 

V-I < 

2.00 (red) 

36/85 

27/72 

2/52 

0/19 

2.00 

< 

V-I 

(reddest) 

30/61 

15/55 

3/56 

0/28 


“ The bulk of the objects within one magnitude below the spectroscopic survey limit are split 
between serendipitously observed objects, which fell within slits designed for brighter targets, and 
objects drawn from the single deeply imaged WFPC chevron number 7. 

* All six objects fainter than 25'* magnitude are serendipitous detections. 


TABLE 3 

Selection Criteria Keys 


Key 

Description 

# 

% 

zsurvey 

redshift sample 

566 

69.2 

serendip 

companions (on shared slitlet) 

121 

14.8 

disk 

elongated galaxies, I < 22.5 

64 

7.8 

align 

mask alignment pointers (mostly stars) 

27 

3.3 

phz 

photometric redshift calibrators 

12 

1.5 

R# 

optical counterparts to known radio source 

9 

1.1 

morph 

moiphologically selected 

7 

0.9 

group A, B 

members of 2 small groups on the sky 

6 

0.7 

hi-z 

dropout candidates 

5 

0.6 

cfrs 

CFRS targets with no redshifts 

1 

0.1 

total 

all spectral targets 

818 

100.0 


The bulk of the spectral targets were selected for the gen¬ 
eral redshift survey program (keyword zsurvey). In addition 
to these primary targets, which were each placed within an 
individual slitlet, we acquired spectra for a set of secondary 
objects (serendip), which also fell within the same slitlets. 
Thus, multiple objects were on occasion observed within a 
single slitlet. This usually occurred by chance, though occa¬ 
sionally by design, where two objects fell together within a 
single default slit length, along or near to the default mask 
angle (40.5°). Some of the serendipitously detected objects 
are extremely faint, and thus we would never have selected 
them as primary candidates for spectroscopic follow-up due 
to the difficulty in getting enough spectral flux into the slitlet. 
These objects (characterized by faint broad band magnitudes 
and very low spectral continuum levels) were typically de¬ 
tected through the fortuitous discovery of an extremely bright, 
isolated emission line. 

A third set of objects {disk) were chosen to allow the acqui¬ 
sition of spatially resolved velocity profiles (Vogt et al. 1996, 
1997,2005). In spite of the misleading key word, these objects 
were not selected for a disklike morphology. They were cho¬ 
sen instead for angular elongation on the sky, equivalent to an 
inclination angle greater than 30° for a late-type spiral disk, 
and for a position angle such that the major axis lay within 30° 
of the Keckn-LRIS mask angle (held at 40.5°). These simple 
criteria allowed the selection of early type galaxies, for exam¬ 


ple {i.e., for which a spatially extended spectrum, composed 
purely of continuum and absorption features but with no emis¬ 
sion was typically observed). Because of the need for high 
S/N spectra, these objects were prioritized only at the bright 
end of the luminosity distribution (Igu < 22.5). 

Every Keckn-LRIS mask which was constructed (see de¬ 
scription below) contained a set of slitlets for science objects, 
and then because the targets were uniformly faint a small 
number of additional boxes (from three to six square slitlets) 
were added for bright, pointlike objects {align) used in the 
mask alignment process. All but two of these objects were 
stars (plus one QSO and one compact galaxy), and as these 
alignment objects were observed along with the science tar¬ 
gets we have spectra for them (albeit at lower resolution, as 
the boxes are two arcseconds wide). 

As mentioned earlier, we prioritized a set of faint, extremely 
blue targets {phz) because we lacked spectra to match HST- 
resolution images of such targets in sufficient numbers to cal¬ 
ibrate photometric redshift models. 

The remaining four subsamples each contain less than 10 
members. For the first set of spectroscopy observations of 
the leftmost one-third of the strip (see Figure |2j, we overlaid 
the detections found in the field at radio wavelengths within 
chevrons 6 or 8-10 (Fomalont et al. 1991) and placed slitlets 
upon those for which an optical counterpart brighter than 24'* 
magnitude existed {R#). An additional set of candidates was 
chosen, particularly in the deeply imaged chevron 7, by eye 
for morphological peculiarities {morph). They each have ei¬ 
ther a diffuse light distribution, or appear to be a potential 
merger in progress, or are a gravitational lens candidate (one 
object). Six objects well within the magnitude limits of the 
redshift sample were included because they appeared to make 
up two small groups in projection on the sky {group A, B). A 
final two objects were included because they fell within the 
region of the strip which had been observed as a part of the 
CFRS {cfrs) redshift survey (Filly et al. 1995), and redshifts 
had not been determined from the CFHT spectra. 

Combining all objects placed within slitlets and the objects 
placed within square boxes for alignment purposes, spectra 
were obtained for a total of 818 unique objects. 

Due to the additional selection criteria, the total set of 
spectroscopic targets is not a strictly random subsample of 
the photometric catalog as defined by our magnitude limits. 
In practice, however, the observed objects sample the range 
of apparent color and magnitude in the photometric catalog 
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fairly evenly in the range (V 606 + f8i4)/2 < 24, as shown in 
Figures 0 through for key structural parameters. Note 
that Simard et al. 2002 and C. N. A. Willmer et al. (2005, 
in preparation) also discuss the detailed selection map in ap¬ 
parent color-magnitude bins. 

8 . MODELING SAMPLE SELECTION 

A set of Monte Carlo simulations was run, selecting objects 
from the imaging catalog purely by aperture magnitudes and 
colors in the fifteen bins to match the fractions listed in Ta¬ 
ble El We began by trimming the spectroscopic sample from 
818 to 782 objects. We first removed six faint, serendipitously 
observed objects which were fainter than the sample limit of 
25 magnitudes. We then removed 30 objects for which Core 
magnitudes could not be measured, due to close proximity 
to WFPC2 chip edges or interfering stellar diffraction spikes. 
The imaging catalog was cut in a similar fashion to 6041 ob¬ 
jects with both magnitudes brighter than 25 and a valid entry 
in the GIM2D catalog. (Note that 947 objects brighter than 
25 were fit by FOCAS but not by GIM2D, which requires a 
larger background area around the object on the sky.) 

Each “mock spectroscopic catalog” contains 782 objects, 
distributed identically to the spectroscopic sample within the 
five color bands and three magnitude bands used in sample 
creation. Figure El shows the distribution of magnitudes and 
colors for the spectroscopic sample and ten mock samples, 
finding reasonable agreement given the scatter in the mock 
catalog distributions. The number of objects within each 
broad magnitude bin is identical for the Core magnitudes and 
colors, by design. 

We observe that the spectroscopic sample contains a slight 
excess of the brightest objects within the 2y‘‘ magnitude band 
(and thus a dearth of the faintest objects), due in part to the 
preferential inclusion of bright stars used in the mask align¬ 
ment process. Within the 24'^’ and 25'* magnitude bands, we 
again see a slight bias towards including brighter objects at 
the expense of fainter ones {i.e., more objects between 23 and 
23.5, or 24 and 24.5, than expected, coupled to less objects 
between 23.5 and 24, or 24.5 and 25). This effect is not larger 
than the variance in magnitudes as measured by different tech¬ 
niques, and becomes washed out in the right hand panel for 
magnitudes derived with the GIM2D package. 

Figure 13 shows the distribution of half light radii, bulge 
fractions, asymmetry indices (also defined as rai\ cf, 
Schade et al. 1995; Simard et al. 2002), a measure of the resid¬ 
ual flux remaining for each object after subtracting a doppel- 
ganger formed by rotating the object image through 180° on 
the sky, and position angle on the sky for the spectroscopic 
sample and ten mock samples. There is a slight trend towards 
including larger objects (by angular size) in the spectroscopic 
survey. The one solid difference between the spectroscopic 
sample and the mock samples, however, is found when ex¬ 
amining the distribution of position angles on the sky. This 
distribution should be relatively flat, though the variance in 
the mock sample distribution indicates the level of scatter. As 
expected, we observe an excess of objects within the spectral 
sample at low position angles. The cause for this trend is the 
prioritized selection of elongated objects (disk classification) 
for which the position angles lie within 30° of the angle on the 
sky at which the Keck-tLRIS masks were placed, causing an 
excess of objects within 10° and 70°. There is a 12% increase 
across this range, a 4tT variation given the scatter in the mock 
catalog distributions. 

Note that this selection effect affects the position angle dis- 



Fig. 7.— The distribution of Core and GIM2D magnitudes and colors, 
for the spectroscopic sample (the hatched histogram, outlined with a heavy 
line) and ten simulated samples drawn from the imaging catalog with the 
magnitude and color selection criteria for the spectroscopic survey but us¬ 
ing no additional criteria (e.g., morphology). There is reasonable agreement 
between the simulated and the actual samples, though a slight trend exists to¬ 
wards brighter objects within each broad magnitude bin for the spectroscopic 
sample Core magnitudes (caused in part, in the brightest magnitude bin, by 
the inclusion of bright stars used as mask alignment objects. 






Fig. 8. — The distribution of angular half-light radii, bulge fractions, Ra 
asymmetry indices, and position angles on the sky, for the spectroscopic sam¬ 
ple (the hatched histogram, outlined with a heavy line) and ten simulated 
samples drawn from the imaging catalog with the magnitude and color selec¬ 
tion criteria for the spectroscopic survey but using no additional criteria (e.g., 
morphology). There is reasonable agreement between the simulated and the 
actual samples for the first three panels. The 4c7 bias in position angle towards 
low values (the dotted horizontal line represents a perfectly even distribution) 
is caused by the preferential selection of elongated objects for which the po¬ 
sition angle of the major axis lay within 30° of the Keck+LRIS mask angle 
of 40.5° (the region marked by dotted vertical lines) within the spectroscopic 
sample. 
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tribution but does not translate into a significantly larger than 
expected fraction of, for example, objects with small bulge 
fractions, as shown in the top panel. Accounting for the 27 
mask alignment objects (e.g., excess stars, which tend to be 
fit with a bulge profile of zero) within the spectral sample, no 
significant differences remain in bulge distributions. In sum¬ 
mary, the selection of a small subset of objects based on addi¬ 
tional factors (e.g., apparent morphology) in the actual spec¬ 
tral sample does not appear to have significantly affected the 
distribution of structural parameters other than position angle 
on the sky. 

We examine the distribution of objects within the complete 
imaging catalog in Figure |3 plotting half light radii, bulge 
fractions, and Ra asymmetry indices against magnitude. In 
order to characterize the distributions, the data were binned 
into a 20 X 20 grid within each panel. Contours were then 
fit to each grid image, running from 5% to 75% of the peak 
values. 

Figure cni contains the subset of data forming the spectral 
sample, plotted in similar fashion. However, the contours in 
Figure are not fit to the spectral sample. Instead, the grid 
images for the total imaging catalog from Figure|9]were con¬ 
volved with a spectral luminosity selection function, and used 
to form a set of contours which reflected the magnitude se¬ 
lection of the spectral sample. This selection function was 
formed by sampling the ratio of counts within the spectral 
sample to counts within the parent imaging catalog, for each 
of the 20 bins formed by the columns of the contour image 
grid (analogous to dividing the counts in the first panel of Fig- 
ure[^by those in the first panel of Figurel^. The effect is to 
multiply each grid point by the fraction of objects selected for 
spectroscopic follow-up within the associated luminosity bin. 

The resultant contours match the distribution of the spectra 
sample quite well, and are in fact indistinguishable from con¬ 
tours drawn directly from the spectral sample data. Figure [m 
shows this comparison explicitly, overlaying the two sets of 
contours against each other. The differences between the dis¬ 
tributions within each panel are not statistically significant. 

9. MASK DESIGN 

The DEEPl spectra were taken as part of a multi-year 
spring season spectroscopy program at the Keck Observatory, 
using the LRIS spectrograph (Oke et al. 1995). Observations 
of the GSS were conducted over nine observing runs on 25 
nights within a five-year period, in conjunction with several 
other spring fields, and thus the process and technique for ob¬ 
ject selection was adapted with passing years as the technol¬ 
ogy of manipulating large field images matured (e.g., memory 
allocation). A total of 36 slit masks was constructed for the 
strip, each containing between 30 and 55 objects. 

An initial set of fifteen LRIS slit masks was constructed in 
1995. We had no equipment which would allow us to create 
slitlets at an arbitrary mask angle at this time, so spectra were 
taken for all objects at the same position angle on the sky 
(40.5°) regardless of the angle of any spatial elongation of the 
object. Individual slitlets were in width and 12" in length, 
and each mask contained roughly 30 of them. The bulk of the 
objects was selected for the general redshift survey, and a set 
of alignment objects was included on every mask. 

The Groth Strip is roughly four arcminutes wide, and so 
we attempted to place objects from only the upper or the 
lower half within a single mask, so as to keep the shift in 
observed spectral range as small as possible from object to 
object. Once objects were selected for the spectral program. 



Fig. 9.— The distribution of (a) number counts, (b) half-light radii, (c) 
Ra. asymmetry indices, and (d) bulge fractions as a function of magnitude 
for all objects within the imaging catalogs, as determined with the GIM2D 
structural analysis package. Dotted lines on every panel separate the points 
into the three bins used in the spectroscopic target selection, with a fourth 
bin for fainter objects extending down to (Veoe + hvA)!'^ - 26. The twelve 
contours overlaid on the last three panels extend from 5% to 75% of peak 
values. 


they were placed upon either one, two, or four masks, depend¬ 
ing upon the faintness of the target. Objects brighter than 23 
in (Veoe + were placed initially on a single mask, those 

between 23'^'^ and 24'^’ magnitude were placed on two masks, 
and the faintest objects, below 24'^‘ magnitude (drawn primar¬ 
ily from the deeply imaged chevron 7), were placed on four 
masks to maximize the amount of spectral flux that we would 
obtain. 

Throughout the spectroscopy program, each mask was 
observed in turn with a blue wavelength grating (typi¬ 
cally 900 1 mm“V5500A) and a red wavelength grating 
(6001 mm“'/7500A), to create a combined spectrum covering 
the range 5000A to 8200A. We planned for 2 x 1500 seconds 
exposure per grating per mask. 

Erom 1996 onward, the minimum acceptable slit-length 
was allowed to drop to 8", as one could extract usable spectra 
reliably with background regions on each side of the object 
of smaller spatial extent due to the overall brightness of the 
continuum level and strong key emission and absorptions fea¬ 
tures. A small fraction of all slits were tilted at an angle (be¬ 
tween 10.5° and 70.5° eastwards of north), in order to trace 
along the major axis of a spatially elongated object or to cap¬ 
ture two objects within one slit. As a result of these changes, 
the average number of objects per mask rose to 45. 

The effect of the distribution on masks along and across 
the entire Groth Strip can be observed in the lower panels of 
Eigure|2l The large, broad peak at -20 arcminutes in the left 
hand panel reflects our initial focus in a region surrounding 
the deeply imaged chevron 7. Note that the spectral fraction, 
the fraction of the imaging catalog objects which were ob¬ 
served with Keck-tLRIS, occasionally rises above a level of 
one, as all spectral targets were counted but the normalization 
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Fig. 10.— The distribution of (a) number counts, (b) half-light radii, (c) 
Ra asymmetry indices, and (d) bulge fractions as a function of magnitude 
for the spectroscopic sample. Dotted lines on every panel separate the points 
into the three bins used in the spectroscopic target selection, with a fourth 
bin for fainter objects extending down to (Veoe + = 26. The spectral 

sample is biased toward the bright end of the magnitude selection bins; the 
fraction of imaging targets contained within the spectral sample drops off sig¬ 
nificantly below V 606 + I 814 = 24, in the two faintest bins (due to the increase 
in the number of possible targets and the decrease in targeting priority below 
this limit). Note that a comparison of bulge fraction histograms between the 
imaging and spectroscopic samples, showing good agreement, can be found 
by comparing the two lower panels in Figure 0 The twelve contours over¬ 
laid on the last three panels were created by convolving the distribution of 
the parent imaging sample, shown in Figure with the fraction of objects 
selected as a pure function of luminosity - they are not fit to these data. Note 
the agreement between the distribution of the spectral sample and the con¬ 
tours, indicating that our complex selection function can be treated as one of 
magnitude without bias in these key parameters. 


was to the limits of the general redshift survey ((Veob + hi 4)/2 
< 24). The right hand panel shows a slight increase in spec¬ 
tral fraction along the top half of the chevrons (chips 2 and 3), 
caused by the distribution of masks along either the top of the 
bottom of the strip. 

10. POSTAGE STAMPS 

A catalog of annotated Veoe and /g^ images has been cre¬ 
ated for the entire spectral sample observed with Keck-i-LRIS 
(Figure El is a representative sample). Each image is over¬ 
laid with the position of the LRIS slitlet at which it was 
observed, and labeled with a set of spectral- and image- 
derived parameters. A single representative panel is shown 
here, and the entire set is available at the DEEP web site at 
http://deep.ucolick.org 


11. SUMMARY 

This paper describes the framework and observational ba¬ 
sis for first five-year phase (DEEPl) of the Deep Extragalac- 
tic Exploratory Probe (DEEP), a longterm, multifaceted study 
of the formation and evolution of galaxies in the distant Uni¬ 
verse. We have obtained complete HST imaging in /814 and 
Veoe of the Groth Survey Strip, a 127 square arcminute region 
in the northern hemisphere sky, and conducted Keck-nLRlS 
spectroscopy of 818 objects in a redshift survey which ex¬ 
tends down to magnitudes of (V 606 + 78i4)/2 = 25. 

Though our selection criteria for the spectroscopic subsam¬ 
ple are complex, and have varied as the survey progressed, the 
bulk of the objects were selected to match the range in colors 
and magnitudes at the bright end {i.e., (V 606 + 78 i 4)/2 < 25) of 
the imaging catalog. We model the complete spectroscopic 
subsample, and show that it can be reproduced statistically 
by making straightforward cuts in magnitude and color in the 
parent imaging catalog. The sole exception to this statement 
lies in the distribution of position angles on the sky, where we 
observe an expected bias towards objects with major axes ori¬ 
ented along the position angle at which the Keck-nLRIS masks 
were aligned. 
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APPENDIX 

APPENDIX A; PARTICULARS OF DEEPl GSS MASK DESIGN 

The DEEPl Groth Strip spectroscopic survey was a multi-year pilot program, designed to find the observational limits of the 
Keck-nLRIS equipment and to explore a variety of scientific programs. For these reasons, as well as the rapid parallel advances 
in computer technologies, the object selection process evolved with time. We present here a detailed overview of the creation 
process for the entire set of 36 masks. 

Sample selection for the first 15 masks was conducted by using a printed set of black and white /gn images of the target 
chevrons, mosaicked by hand into a single, extended map. Transparent overlays were coded to indicate object Core magnitudes 
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Fig. 11.— A set of contours fit to the distribution of the spectral sample (solid lines) are compared to those from Figure fTol (dotted lines), which were 
creating by convolving a set of contours fit to the distribution of the total imaging catalog with the fraction of objects observed in the spectral sample as a 
function of magnitude. In the first panel, note that the small dotted contour in the region of high values falls at the lowest contour (5% of peak values), and 
is not statistically significant. Both distributions of bulge fractions (second panel) show a strong peak for objects with extremely low bulge fractions, while the 
distributions of half light radii (third panel) are broader. There is good agreement between the two fits in every panel. 



Fig. 12.— A representative set of 1^06 and HST+WFPC2 images of objects within the spectral sample. The solid lines on each postage stamp image 
indicate the position of the KECK+LRIS slitlet upon the object. Each object is labeled by name, followed by total /gi 4 GIM2D magnitude, aperture /gi 4 Core 
magnitude, absolute B-band magnitude, and bulge fraction; the lower line contains redshift, V 606 “f 8 i 4 GIM2D color, aperture V 606“^814 Core color, half light 
radius in arcseconds, and asymmetry parameter Ra- 


and colors. This first set of masks was placed to overlap with the deeply imaged chevron 7 at the northeast end of the strip (fourth 
chevron from the left end, in Figure|2j, and extended from chevron 6 down to chevron 10. The first group of four masks (GSS1 A, 
GSSIB, GSS 1C, and GSS ID) had a single position on the sky and contained objects selected from chips 2 and 3 of the targeted 
chevrons, while the second and third sets of four (GSS2-ABCD, GSS3-ABCD) contained objects selected from chips 3 and 4 and 
were offset toward the lower edge of the Groth Strip. The final three masks (GSSB-abc) were back-up program masks designed 
for use in poor weather (low transparency conditions), with 1.2'^ width slitlets. Each mask contained the brightest candidates 
from the GSSl, GSS2, and GSS3 series respectively, and were filled with additional bright objects in the range (Veoe + huV'^ 
<23. 

A few serendipitous spectra would be discovered during the spectral reduction process, and the entire set of radio and group 
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objects was selected at this time, as well as several morphological candidates. Note that no objects were selected according to 
the disk, cfrs, or phz key words for these early mask sequences (though the bluest and reddest objects were prioritized within the 
general redshift survey). 

Beginning in 1996, masks were designed using interactive, computerized selection to select objects and to place them on the 
masks (NPV’s disktool IRAF package). Following target and alignment star selection, all masks were modeled with ACP’s 
ucsclris IRAF package to account for specific observational constraints {e.g., anamorphic corrections) and to create machine 
(CAD) instructions for milling the mask locally in the UCSC instrumentation shops. The pattern of placing objects on one, two, 
or four masks according to brightness was maintained. In addition, any objects which were re-observed because of insufficient 
initial flux levels {i.e., no redshift could be assigned) were placed on as many masks as were available, and objects with the 
reddest colors were automatically placed on at least two masks regardless of brightness (spectral reductions proved to be difficult 
for these objects, due to the presence of very little other than wide, diffuse absorption features). After the first year a small 
fraction of all slits were tilted at an angle {0 < 30°) relative to 40.5°, in order to trace along the major axis of a spatially elongated 
object or to capture two objects within one slit. For the disk candidates, exposure time was doubled from one mask to two to 
increase S/N levels for rotation curve extraction. 

In 1996, the thirteen unused masks were updated with several elongated disk candidates, several additional morph candidates, 
and a set of the phz objects with extremely blue colors. Because KPNO UBRI photometry had been obtained for the northeast 
portion of the strip in the intervening year, these data were used as well in the relative prioritization of the 17 phz candidates. 
The masks were updated in place, with original names retained. (Only masks GSSl-A and GSSl-B were used in 1995, due to 
bad weather.) In 1997 we created a new set of five masks (GSS4-ABCDE), shifting our default position along the top of the 
strip slightly to the right to include chevron 11. We ceased placing general targets on the masks which were fainter than (Veoe 
+ I 8 i 4 = 24, as we found them to demonstrate unlikely success for redshift identification. Five UB dropout candidates (with 
photometric redshifts z > 2.5) were identified from the KPNO UBRI data and placed on four masks each. In 1998 seventeen 
additional masks were created (GSS5-ABC, GSS6-ABCDEF, and GSS7-ABCDEEGH). We began to employ tilted slits not only 
to trace the major axes of individual spatially extended objects but also to place secondary, and occasionally tertiary, serendip 
objects which were offset from primary targets by 0 < 30° from the mask angle, to increase efficiency. 

The GSS5 series was designed to re-observe objects from within chevrons 6 through 11 for which we still did not have redshifts, 
in spite of repeated observations. We prioritized objects in the following order: (1) previously observed targets with Ig^ < 23 
and no secure redshift, (2) previously observed sources with Ig^ < 24 and no secure redshift, (3) previously observed sources 
with Igi 4 > 25 and no secure redshift, (4) high redshift (U and B dropout) candidates determined from ground based UBRI data 
or from existing GSS spectra which showed promise (i.e., no features indicating lower redshifts), and (5) CERS survey objects 
without redshifts. The GSS6 series was placed at the southwest end of the strip. Two new pointings were established, and at each 
one two normal masks (ABDE) and one back-up mask (CE) with bright objects was created. The disk, phz and extremely red 
objects were prioritized in the sample selection process. The GSS7 series was placed in the middle region of the strip, designed 
to sample the bright end of the luminosity function and bridge the spatial gap between the northeast and southwest sampled 
portions. Eour new pointings were created, and a normal mask (ACEG) and a back-up mask (BDEH) with bright objects were 
made for each one. The disk, phz and extremely red objects were again prioritized in the sample selection process; in addition as 
there would be only one mask exposure per position under optimum conditions, we limited the general redshift survey to (Veoe + 
l8i4)/2 < 23. 

In 1999 ten additional masks were created (GSS6-GHIJ and GSS7-IJKLMN). They were placed in sets of two at five of the 
six positions established in 1998 (the sixth pointing, home of masks GSS7-EE was not doubled because it had not been used in 
1998). 

We note two factors of interest to anyone studying interactions, pairs of objects, or large scale structure: (a) the repeated 
placement of fainter objects on a series of masks all centered at the same location, conducted to increase the total spectral flux 
and thus the probability of redshift determination for such objects, and (b) the requirement of a few (four to six) arcseconds 
of slit being available to each side of each object for background light removal. These both affect the object-to-object angular 
distribution within the spectral sample on small size scales (9 ~ 10")- The first effect will further result in a sparser sampling of 
near neighbors on the sky for the fainter objects within the spectral sample. 

APPENDIX B: CONVERSION EACTORS BETWEEN AB AND VEGA MAGNITUDES 

The convolution between filter responses and galaxy SEDs followed Eukugita et al. (1995), by resampling filters and spectra 
to the same dispersion (1 A), using parabolic and linear interpolations respectively. Eor the HST /gi 4 and Veoe filters, efficiency 
curves for WEPC2 were downloaded from the Space Telescope Science Institute web site and curves for the different CCDs were 
averaged. The Vega spectrum used to calculate the conversion to AB magnitudes is the model atmosphere calculated by Kurucz, 
distributed with the Bruzual & Chariot (2003) galaxy evolution synthesis package. The AB spectrum is simply a flat spectrum in 
E(i/) converted into wavelength space (e.g., Eukugita et al. 1995). 
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TABLEB4 

Conversion Factors between AB and Vega Magnitudes 


Transformation Between AB and Vega Systems 

V606ab 

= V606\/e^fl + 0.096 

m^AB 

= /814v,sa +0.417 

(V606-/814)ab 

= (V606-/814)v<.^„-0.321 
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